m 




Spectroscopy with dynamical Chirally Improved 
quarks 



Christof Gattringer", C. B. Lang", Markus Limmer ', Thilo Maurer^, Daniel Mohler* 
and Andreas Schafer'' 

" Institutfilr Physik, FB Theoretische Physik, Universitdt Graz, A-8010 Graz, Austria 
* Fakultdtfiir Physik, Universitdt Regensburg, D-93040 Regensburg, Germany 



E-mail: christof . gattringer@uni-graz . at[ christian . lang@uni-graz . at 


markus . limmer@uni-graz . at 


thilo . maurergphysik . uni -regensburg . de 


daniel . mohler@uni-graz . at 






andreas . schaef er@physik . uni -regensburg . de 





We present recent results of our dynamical simulations with Chirally Improved fermions and 
report on new developments in the determination of excited light-quark meson states using inter- 
polators constructed by applying covariant derivatives on Jacobi-smeared quark sources within 
the framework of the variational method. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Spectroscopy with CI quarks 



Daniel Mohler 



1. Introduction 

Simulations with full chiral symmetry using Overlap fermions are still in their infancy. Cur- 
rently the computational cost of such simulations prohibits large-scale simulations of fermions 
fulfilling the Ginsparg-Wilson (GW) relation exactly. Therefore, we are working with an approxi- 
mate algebraic solution to the GW relation [|l|,|^. These so-called Chirally Improved (CI) fermions 
are still computationally demanding as the fermion action includes several hundred terms ranging 
up to three links distances. 

Simulations with CI fermions in the quenched approximation have demonstrated good chiral 
and scaling behavior [Q]. Interpolators constructed from Jacobi-smeared sources [Q, ^] of different 
width as well as interpolators containing derivatives have been used to extract the spectrum of light 
mesons [|6[ ^] within this approximation. For the extraction of excited states the variational method 



[g, ^] has been used, which also enables one to separate contributions from ghost states []10|]. 

Here we present a progress report of ongoing dynamical simulations with CI fermions. The 
focus will be on the mass spectrum and we will present preliminary results for both mesons and 
baryons. Preliminary results on the simulation and hadrons masses have previously been presented 
in [11, O] . For a discussion of other recent results for excited states see [ 1 3 ] . 



2. Details of the simulation 

For our simulations with dynamical quarks we use the Liischer-Weisz gauge action and two 
mass degenerate CI fermions. We use a a standard HMC algorithm with mass preconditioning [[I^ 
and a mixed-precision inverter. Our action also incorporates one level of stout smearing. Table |l] 
shows the simulation parameters for our configuration ensembles A-C (lattice size 16^ x 32). 



ensemble 




mo 


HMC time 


^[fm] 


m;t[MeV] 


A 


4.70 


-0.050 


591 


0.1507(17) 


526(7) 


B 


4.65 


-0.060 


1108 


0.1500(11) 


469(4) 


C 


4.58 


-0.077 


1046 


0.1440(11) 


318(5) 



Table 1: Run parameters, lattice spacing and pion masses for ensembles A-C of lattice size 16^ X 32. The 
scale has been set using the Sommer parameter ro.exp = 0.48 fm. 



The lattice spacing has been determined from the static quark potential (for hypercubic smeared 
configurations). We fit to the potential 

VL{r)=A+- + or + C3AV{r) , (2.1) 

where AV{r) = [7] — 7 contains the perturbative lattice Coulomb potential [i] . Using a value of 
''o.exp = 0.48 fm for the Sommer parameter we can then determine the lattice spacing as follows: 

Figure |l] shows example fits to the logarithm of the Wilson loop expectation value \nW{r,t) and to 
the potential for ensemble C. 
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t r 

Figure 1: L.h.s.: Fit to lnW{r,t) in the range 4 < f < 7. R.h.s: Fit to the potential in the range 1 < r < 7. 
The data points have been connected by straight hnes to guide the eye. 



3. Quark sources for the variational method 

3.1 Variational method 

For the extraction of hadron masses we use the variational method [§, One uses a matrix 
of correlators projected to fixed (here: zero) spatial momentum 

C{t)ij = "L{^\Oi\n){n\o]\0) , (3.1) 

n 

and solves the generalized eigenvalue problem 

C{t)vk = h{t)Cito)n, A,(0-e-'^*(l + ^(e-'^*)) . (3.2) 

At sufficiently large time separation each eigenvalue receives contributions from only a single mass. 
At the same time the eigenvectors can serve as a fingerprint to identify the states when followed 
over several f-values. For a recent discussion of the generalized eigenvalue problem in that context 
see [0]. 

3.2 Quark sources 

We construct interpolators from Jacobi-smeared ^ quark sources S of different width, with 
the same choice of the parameters K and A'^ as given in |^. 

N 

S = MSq with M=Y,T^H" and (3.3) 

n=0 

H{n,m) = £ (uj{n,0)5{n + j\m)+Uj{n-j,0y 5{n-j,m)) . 

denotes a point source. Combinations of such approximately Gaussian sources allow for nodes 
in the interpolating operators while using fewer quark propagators than in other approaches. In 
addition to these sources we also use derivative quark sources W^/, by applying covariant derivatives 
to the Gaussian smeared sources: 

Wrf, = DiS , Di{x,y) = Uiix,0)5ix + i,y) - Ui{x-i,0)'^ 5{x-i,y) . (3.4) 
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Figure 2: 1st and 2nd excitation of n. The shaded region indicates the systematic uncertainty due to the 
choices of fit-range and interpolator combination. 



4. Spectrum results 

In this section we present results from quenched lattices as well as preliminary results from 
the simulation with dynamical CI quarks. 

4.1 Results from quenched lattices 

Figure § shows results for the pion channel based on 99 quenched configurations. A clear 
signal can be obtained for both the first and the second excitation. The results from the set of inter- 
polators containing both Gaussian and derivative sources enable us to identify the second excited 
state which (on the same dataset) cannot be observed with the Gaussian sources alone. For more 
details, some further plots and results from other channels see [^. 

4.2 Results from dynamical lattices 

For the dynamical CI lattices, the lattice size is 16^ x 32, corresponding to a spatial extent of 
(2.4 fm)^ X 4.8 fm. We analyze every fifth configuration and shift the source positions for consecu- 
tive configurations to reduce autocorrelation. All data presented are of a preliminary nature. Table 
^provides an overview of the current ensembles. 

4.2.1 Results for baryons 

Figure ^ shows results for the mass of the (positive parity) nucleon ground state. The empty 
symbols indicate a nucleon made from valence quarks that are heavier than the sea quarks while the 



ensemble 


# conf.s 


a [fm] 


wjawi [MeV] 


mjt [MeV] 


amji L 


A 


100/100 


0.1507(17) 


43.0(4) 


526(7) 


6.4 


B 


100/200 


0.1500(11) 


35.1(2) 


469(4) 


5.8 


C 


100/200 


0.1440(11) 


15.0(4) 


318(5) 


3.8 



Table 2: Dynamical CI runs. # conf.s refers to the number of (independent) configurations analyzed for 
mesons or baryons respectively. Values in physical units have been obtained as described in Section 2. 
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Figure 3: A: Mass of the (positive parity) nucleon ground state; B: Mass of the negative parity ground state 
and 1st excited state. 
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Figure 4: A: Ground state in the Delta channel in physical units; B: Same plotted in units of the nucleon 
mass. 



full symbols indicate the fully dynamical, unitary point for the respective runs. The error bars are 
purely statistical, also including the statistical error in the determination of the scale. Systematic 
errors from varying fit ranges and combinations of interpolators are not included in the error bars. 
Since we have set the scale individually for each of the parameter sets this amounts to a mass 
dependent renormalization scheme. We therefore cannot use ChPT extrapolations valid in the mass 
independent scheme. More runs at different masses and volumes will be necessary for reliable 
extrapolations to the physical points. 

Figure |3|B shows the ground and first excited state masses observed in the negative parity 
channel. While the statistical errors for the different runs show a deviation a more careful analysis 
of the systematics due to the choice of interpolators and fit ranges will be necessary. The quality 
of the data for run C was not sufficient to fit these states at the unitary point, therefore it has been 
omitted from the plot. 

Figure ^ shows the ground state mass in the A-channel plotted in physical units and in units of 
the nucleon mass. The A++ is a resonance, which complicates the interpretation of the ground state 
signal. To get a clear interpretation an analysis of this state in different volumes will have to be 
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performed. We observe that when plotting the A++ in units of the nucleon mass the finite volume 
effects at small pion masses cancel, thus leading to consistent data for all three runs. 
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Figure 5: A: Ground state in the 1 channel (p-meson); B: Ground state signal in the 0^^ channel for 
interpolators 9, 10 in the notation of [^. 



4.2.2 Results for mesons 

For the mesons the set of interpolators containing derivative sources described in detail in 
has been used. As can be seen in Table |2[ only a subset of configurations was used for this 
preliminary analysis. Figure ^ shows some results for the ground state of the p meson. Figure |5|B 
displays results for the lowest energy level in the isovector scalar ao channel for one possible set 
of interpolators. We point out, however, that our data show a large operator dependent variation. 
Most quenched results' indicate that the ground state in this channel might be the ao(1450), in 
agreement with a tetraquark interpretation of the a()(980). Results from dynamical simulations 
lead to somewhat lower states, possibly due to the presence of nrj scattering states. A systematic 
analysis will have to take into account several volumes to discriminate scattering states from bound 
states. One possible interpretation of our results would be a crossing of the lowest energy level at 
the quark mass where the nrj state becomes lower than the ao(980). 



5. Conclusions 

We presented a progress report on spectroscopy with dynamical CI fermions. We demonstrated 
that good signals can be obtained for (most) ground state mesons and baryons. Emphasis has 
been put on describing the methods and we are currently refining these by experimenting with 
different gauge link smearings. It would be desirable to perform further runs with the CI action 
to control both discretization and finite volume effects and to be able to use a mass independent 
renormalization scheme. 

'Notice that the analysis of the quenched data is made difficuU by artifacts, so called "ghosts". 
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